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■ We study the effective monopole action obtained in the Maximal AbeHan projec- 
• tion of quenched SU(2) lattice QCD. We determine the quadratic part of the lattice 
, action using analytical blocking from continuum dual superconductor model to the 
\ lattice model. The leading contribution to the quadratic action depends explicitly 

^ I on value of the monopole condensate. We show that the analytical monopole action 

■ matches the numerically obtained action in quenched SU(2) QCD with a good ac- 
I ■ curacy. The comparison of numerical and analytical results gives us the value of the 

CN . monopole condensate in quenched SU(2) QCD, rj = 243(42) MeV. 

^ PACS numbers: 11.15.Ha,12.38.Gc,14.80.Hv 
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O ; I. INTRODUCTION 

o 

■ The dual superconductor mechanism jl[ is one of the most promising mechanisms invented 
^ . to explain the confinement of color in non-Abelian gauge theories. The basic element of 

\ this mechanism is the existence of specific field configurations - called Abelian monopoles - 
I I in the QCD vacuum. The monopoles are identified with the help of the Abelian projection 
method 13], which uses the partial gauge fixing of the SU(N) gauge symmetry up to an 
^ . Abelian subgroup. The Abelian monopoles appear naturally in the Abelian gauge as a 
^ \ result of the compactness of the residual Abelian group. 

Various numerical simulations indicate that the Abelian monopoles may be responsible for 
the confinement of quarks (for a review, see, e.g., Ref. |3]). The Abelian monopoles provide 
a dominant contribution to the tension of the fundamental chromoelectric string ^ ^, 0| . 
In Ref. 01 it was qualitatively shown that the monopole condensate is formed in the low 
temperature (confinement) phase and it disappears in the high temperature (deconfinement) 
phase. The energy profile of the chromoelectric string as well as the field distribution inside 
it can be described with a good accuracy by the dual superconductor model 0, H, |^ 

There were various attem pts to determine the dual lagrangian and the values of its 
couplings millilnlHHlIi- A simplest version of the dual superconductor model for 



X 
J3 



SU(2) gauge theory contains three independent parameters: the mass of the monopole, M$, 
the monopole charge, and the value of the monopole condensate, 77. In Ref. [s*] the SU(2) 
string profile was compared with the classical string solution of the dual superconductor in 
continuum and the mass of the dual gauge boson, Mb = gt], and the monopole mass were 
shown to be equal^, Mb ~ M$ ^ 1.3 GeV. These values are close to the results of other 
groups. 



In this paper we quote the first set of parameters of Ref. ^8] which is self-consistent. 
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The value of the monopole condensate derived from the chromoelectric string analysis 
of Ref. is ?7 = 194(19) MeV. In this paper we determine the value of the monopole 
condensate from the effective monopole action obtained in the numerical simulations of 
quenched SU(2) QCD. Our strategy is the following. We relate the lattice monopole model 
on the lattice with the continuum dual superconductor model using the approach of blocking 
of the continuum variables to the lattice proposed in Ref. Generally, this method allows 
to construct perfect lattice actions and operators in various field theories. In particular, 
this method was used in Ref. \16^ for the quenched SU(2) QCD at high temperatures to 
study the dynamics of the static monopoles. The lattice monopole action obtained with the 
help of such a blocking depends on the parameters of the original continuum model. The 
comparison of the analytical form of the lattice monopole action with the corresponding 
numerical results allows in general to fix the parameters of the continuum model. In this 
paper we concentrate on the determination of the monopole condensate in the quenched 
SU(2) QCD in the Maximal Abelian projection [l7| . 

The plan of the paper is the following. In Section IH] we propose the method of blocking 
from continuum to the lattice of the monopole currents in four dimensional space-time. 
We compute the quadratic part of the monopole action analytically in Section IIIIl while 
the numerical computation is done in Section HVl In Section |3 we compare the numerical 
data with analytically calculated action and fix the value of the monopole condensate. Our 
conclusions are presented in the last Section. 



II. BLOCKING FROM CONTINUUM IN FOUR DIMENSIONS 

The method of blocking of continuum variables to the lattice 3, 3 constructs the lattice 
model (at given finite lattice spacing b) starting from a model in continuum. The essence 
of this method is simple. Consider, for example, the blocking of the topological variables, 
such as the monopole charge in three space-time dimensions In 3D the monopoles 

are instantons characterized by their positions and the magnetic charges. Suppose, that the 
dynamics of these monopole charges in continuum is described by a Coulomb gas model 
with two parameters, the fugacity, (, and the monopole charge, g. Let us superimpose a 
cubic lattice with the lattice spacing on a particular configuration of the monopoles. Each 
of the lattice 3D cells can be characterized by an integer magnetic charge it contains. Thus 
we can relate the continuum configuration of the monopoles to the lattice configuration 
characterized by magnetic charge inside each cell (see Figure ^ for an illustration). Next 
step is to construct a "lattice quantity" (for example, the absolute value of the magnetic 
charge inside a 3D cell) and calculate analytically the average of this quantity over all 
configurations of the continuum monopoles. The value of this averaged quantity would 
depend on the size of the cell, b, and on the parameters of the continuum model {i.e., on 
( and g). Similarly, one can study numerically the same quantity in a purely lattice model 
(i.e., in the dimensionally reduced quenched SU(2) QCD as in Ref. |l6|), and relate both 
numerical and analytical results for the density with each other. Since the averaged density 
depends on the scale b, the fitting of the numerical results to the analytically obtained 
formula gives an information about the parameters of the continuum model, ( and g. The 
fitting also provides an information about the self-consistency of this approach, or, in other 
words, about the validity of the description of the lattice quantities by the continuum model. 

Therefore this method allows us to describe the lattice observables by the continuum 
model. In Ref. 3| the blocking was performed for the monopoles in 3-D which are the 
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(a) (b) 



FIG. 1: Blocking of the continuum monopoles to the lattice in (a) 3-D and (b) 4:D. In 3-D the 
charge corresponding to the lattice cube C is given by the total magnetic charge of the continuum 
monopoles inside this cube. In 4Z) the charge is proportional to the linking number of the monopole 
trajectory, k, with the surface of the three-dimensional cube C. 



instanton-like objects. Below we generalize this approach to the four dimensional case. 

The partition function of the dual superconductor can be described in terms of the 
monopole trajectories as follows: 

Z^on = jvk j VB expj- y d^x [^Flu + ik^{x)B^{x)\ - , (1) 

where F^^ = d^By — di,B^ is the field stress tensor of the dual gauge field B^, and Sintik) is 
the action of the closed monopole currents fc, 

k,{x) = ydr^^5(^)[x-x(r)]. (2) 

Here the 4D vector function x^(r) defines the trajectory of the monopole current. In Eq. 
the integration is carried out over the dual gauge fields and over all possible monopole 
trajectories (the sum over disconnected parts of the monopole trajectories is also implicitly 
assumed) . 

The action in Eq. (Q) contains three parts: the kinetic term for the dual gauge field, the 
interaction of the dual gauge field with the monopole current and the self-interaction of the 
monopole currents. The integration over the monopole trajectories gives the lagrangian of 
the dual Abelian Higgs model [lO |: 



ZfYicyri^ OC Z 



DAHM 



jv^ j VBexp{-J d'x[^F'^, + ^\{^, + ^B,)<^\' + (3) 



where $ is the complex monopole field. The self-interactions of the monopole trajectories 
described by the action Sint in Eq. (jT)) lead to the self-interaction of the monopole field $ 
described by the potential term V{^) in Eq. 
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Now let us embed the hypercubic lattice with the lattice spacing b into the continuum 
space. The three-dimensional cubes are defined as follows: 



Cs,fM = ^b(^Su - ^ j < <b(^s,, + for 1/ 7^ /i ; and = 6s^| 



(4) 



where is the dimensionless lattice coordinate of the lattice cube Cs,fj. and x,^ is the con- 
tinuum coordinate. The direction of the 3D cube in the 4D space is defined by the Lorentz 
index /i. 

As in the three-dimensional example described above, let us consider a configuration of 
the monopole currents superimposed on the lattice (@)). The monopole charge Kc inside 
the lattice cube C^^^ is equal to the total charge of the continuum monopoles, k, which pass 
through this cube. Geometrically, the total monopole corresponds to the linking number 
between the cube C and the monopole trajectories, k (an illustration is presented in Fig. P). 
The mutual orientation of the cube and the monopole trajectory is obviously important. 
The corresponding mathematical expression for the monopole charge Kc inside the cube C 
is a generalization of the Gauss linking number to the four dimensional space-time: 



Kc{k)^E{dC,k) = ljd'^ jd'ye,,^p^f^{x)k^{y)dpV^'\x-y) 

= 4^ / 1 d'y e,.„, S2^(a:) k^{y) ^^-^ . (5) 



Here the function S^^(x) is the two-dimensional 5-function representing the boundary dC 
of the cube C. In general form it can be written as follows: 



= dndr, - x{f)] , (6) 



where the four dimensional vector x{f) parameterizes the position of the two-dimensional 
surface E. The function T>'^^'> in Eq. (0) is the inverse Laplacian in four dimensions, 
SpD^^^x) = 5^^\x). It is obvious that the lattice currents Kg^f^ are closed 

d'K = 0, (7) 

due to the conservation of the continuum monopole charge, d^k^ = 0. In Eq. ((7j) the 
symbol d' denotes the backward derivative on the lattice. We present a proof of Eq. ^ in 
Appendix ^ 

Let us rewrite the dual superconductor model Q in terms of the lattice currents Kc, 
Eq. (0). To this end we insert the unity, 

1= E ll^{Kc-IL{dC,k)), (8) 

Kc&^ C 

into the partition function (here 6 represents the Kronecker symbol). Then we integrate 
the continuum degrees of freedom, k^ and B^, getting the partition function in terms of the 
lattice charges Kc- The simplest way to do so is to represent the product of the Kronecker 
symbols in Eq. (jH)) in terms of the integrals, 

oo 

1= XI [n / ^^^] ^w{^Y.^cKc-^ I d^x fc^(a;)5^(^^; x)} , (9) 

Kc&Z C C 
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where 



c 



(10) 



To derive Eqs. from Eq. (jHl) we used relation ©. 

Substituting Eq. © into Eq. (0) we get: 



l^k I VB [U 

1 



exp 



c 



—Fl^ + ik,{x) ( B,{x) + x] 



Sint{k) 



(11) 



One can see that the substitution of the unity © effectively shifts the gauge field in the 
interaction term with the monopole current, B^ — > i?^ + i?^. Therefore the integration over 
the monopole trajectories, k^, in Eq. is very similar to the integration which relates 
Eq. ((H) and Eq. (jH)). Thus, we get: 



2mon Z^AHM 



P$ / 5^ [J] / dOc 

Kc&^ C 



exp 



c 



Summarizing this Section, we rewrite the continuum dual superconductor model in terms 
of the lattice monopole currents, K: 



Zdahm — > e 



E 



(K) 



(13) 



where the monopole action is defined via the lattice Fourier transformation: 



Smon 

(K) 



VOc exp<^ -S{e)+i{e,K) 



(14) 



and the action of the compact lattice fields 9 is expressed in terms of the dual Abelian Higgs 
model in continuum: 



e-^W= VB exp\- d'^x 



4 2 M^^^2 



d^ + t{B^ + B^i9) $ +Vi^) . (15) 



Eqs. fll()piH|14pi5|l are the main result of this Section. 



III. QUADRATIC PART OF MONOPOLE ACTION 



An exact integration over the monopole, $, and dual gauge gluon, B^, fields in Eq. flT3)) 
is impossible in a general case. However, in this paper we are interested in the quadratic 
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part of the monopole action which is dominated by the contribution of the one dual gluon 
exchange. Therefore we do not consider the effect of the fluctuations of the monopole field $, 
which lead to the higher-point interactions in the effective monopole action^ ^M- Effectively, 
the neglect of the quantum fluctuations of the monopole field corresponds to a mean field 
approximation with respect to this field, $ ($). In this case the AHM action becomes 
quadratic and Eq. (fTHjl can be rewritten as 

e-'^'^= J^Bexp{-Jd'x[-^X^ + ^-^[B, + B,ie)y]}, (16) 

where rj = \ ($) | is the monopole condensate. 

The Gaussian integration over the dual gauge field can be done explicitly. In momentum 
space the effective action (up to an irrelevant additive constant) reads as follows: 



where B^{9,p) is related to the field B^{9,x), given in Eq. (fTUj) . by a continuum Fourier 
transformation: 

B,ie,p) = ^ J2[p^ S^aQaipb) - P,Pa Qa(p&)]e-^'(^'^) 9s,a , (18) 



with 



sinx^/2 



To get Eq. ([T^ from Eq. (|Tnjl we notice that 

\e,,^p^i'l{x) = dy,,V^^{x) , (20) 



where is the characteristic function of the lattice cell Cs^^i- Namely, the characteristic 
function of the 2>D cube with the lattice coordinate and the direction a is 

V^{Cs,a, x) = 5{x^ - hs^) J] Q{h{s^ + 1/2) - x^) ■ Q{x^ - h{s^ - 1/2)) , (21) 

where G(x) is the Heaviside function. The Fourier transform of the function ()21|) is 

l^M(Cx,a,p) = 5m," QM e-^'(^'^) . (22) 

Substituting Eq. (fTHj) into Eq. (fTTj) and changing the momentum variable, q = bp, we get 
the following expression for the quadratic action: 

S{0) = —— ^ ^ 0s,a^~),^^,9s>a0s>,a' , (23) 



^ The fluctuations of the monopole fields and their effect on the blocked monopole action will be considered 
in a subsequent paper. 
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where 

= / ^ ^'^-'-Jf" g,(,)g,,(g) e-(^'-) . (24) 

Here we have introduced the dimensionless parameter 

fx = MbB. (25) 

Next step is to substitute Eq. into Eq. (fT^ and integrate over the variables 9 to get 
the quadratic monopole action: 



We could not find an explicit form for the operator and therefore we calculate it 
in the n —>■ oo limit. This limit corresponds to large values of b which are consistent with 
the quadratic form of the monopole action . The details of the calculation are given in 
Appendix [Bl and the result is: 

+-ga{{s- s')^) + —^Ho^iis- s')^)] , (27) 

cyclic cyclic ^^^'^ 

Here Da\s±) is the three-dimensional Laplacian acting in a timeslice perpendicular to the 
direction a, Ss is the Kronecker symbol, = T)^^\s) is the one-dimensional Laplacian 
operator (double derivative) defined in Eq. ()B15|) . T{a, x) is the incomplete gamma function 
and tuv is an ultraviolet cutoff. In Eq. (j77j) exponentially suppressed corrections of the 
order 0{e~^°'^^^'^) are omitted. 

Inverting the operator (P7j) and expanding it in inverse powers of fi we get the quadratic 
operator S in the monopole action 



(4^ T^-.' T^'o.' T^-o.' - srv-'n^ D^i) + o(/i-^)] . (28) 



where = T>a \ F = r(0, t[/yM^ 6^). The operator expansion in Eq. (PHj) is written in a 
symbolic form. 



IV. MONOPOLE ACTION IN QUENCHED SU(2) QCD 



Having determined the action of the blocked monopoles analytically, we are going to the 
same in the quenched SU(2) QCD using numerical calculations. We simulate the quenched 
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SU(2) gluodynamics with the lattice Wilson action, S{U) = — f ^pTrf/p, where (3 is the 
coupling constant and Up is the SU(2) plaquette constructed from the link fields. All 
our results are obtained in the Maximal Abelian (MA) gauge ^3] which is defined by the 
maximization of the lattice functional 

R = Y,^^{^3U{s,fi)a3U\s,fi)') , (29) 

with respect to the SU(2) gauge transformations U{s,fi) — > U{s,fi) = Q{s)U {s, fi)Q'^ {s + fi) . 
The local condition of maximization can be written in the continuum limit as the differential 
equation {df^+igA'^^){Aj^—iA'j^) = 0. Both this condition and the functional (j29p are invariant 
under residual U(l) gauge transformations, ^l^^'^^{to) = diag(e*'^'''^-', e"*'^'-'*^). 

After the gauge fixing is done we get the Abelian variables applying the Abelian projection 
to the non-Abelian link variables. The Abelian gauge field is extracted from the SU{2) link 
variables as follows: 

where u{s,^) = exp(i6'(s, /i)) represents the Abelian link field and c(s, /x) corresponds to 
the charged (off-diagonal) matter fields. The Abelian field strength 9^u{s) G (— 47r, 47r) 
is defined on the lattice plaquettes by the Abelian link angle 9{s,fi) G [— 7r,7r) as follows: 

V(s) = o{s, n) + e{s + fi,u)- e{s + i>,fi)- e{s, u). 

To construct the Abelian monopoles we decompose the field strength d^v{s) into two 
parts, 

d^,{s) = d^,{s) + 2'Km^,{s), (31) 

where 9^ij{s) G [— vr, vr) is interpreted as the electromagnetic fiux through the plaquette 
and m^^{s) can be regarded as a number of the Dirac strings piercing the plaquette. The 
elementary (i.e., defined on the 1'^ lattice cubes) monopole currents are determined by the 
DeGrand-Toussaintjl^ formula: 

kf,{s) = ^e^^p^d^mp^{s + fi), (32) 

where d is the forward lattice derivative. The elementary monopole current is defined on 
a link of the dual lattice and takes values 0, ±1, ±2. Moreover the elementary monopole 
current satisfies the conservation condition by construction, 

d'M') = ' (33) 

where d' is the backward derivative on the dual lattice. 

Besides the elementary monopoles one can also study the so called extended monopo- 
les jiO]. The extended monopoles are usually used to define the monopole current on a cube 
of a large size without getting artificial lattice corrections of the order of the lattice spacing, 
a. The extended monopole is defined on a sublattice with the lattice spacing b = na. The 
explicit construction of the extended monopoles corresponds to a block spin transformation 
of the monopole currents with the scale factor n, 

n-l 

k^-^\s) = J2 + {n- l)/i + ii> + JP + la) . (34) 

i,j,l=0 
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The spacing a of the original lattice and, consequently, the artificial lattice corrections (which 
are of the order of 0(a)) can be arbitrarily small while the size of the blocked monopole can 
be fixed in physical units. In our studied we have studied n = 2, 3, 4, 6, 8 blocked monopoles 
on 48^ lattice. 

Applying consecutively the gauge fixing, the Abelian projection and using formula ()34|1 
one can construct the Abelian monopole ensemble for any ensemble of the non-Abelian 
fields of quenched SU(2) QCD. Then using an inverse Monte-Carlo method one can get the 



effective monopole action. The details of this procedure can be found in Refs. |18l . |21L |2^ . 
In our simulations we have used 200 configurations on 48^ lattice. The Maximal Abelian 
gauge was fixed with the help of the standard iterative procedure. 

In general, the monopole action, 5^°°, can be represented as a sum of the n-point {n > 2) 
operators Si, Ref. ||18i. .21i]: 



S[k] = Y,g,S,[k] 



(35) 



where fi are coupling constants. In this paper we adopt only the two-point interactions of 
the form Si ~ k^{s)k^i{s') which works well at large values of h. Using the inverse Monte- 
Carlo method we calculate the monopole action parameterized by 27 couplings gi. The 
maximal distance between the interacting currents in this action is 3 in units of the blocked 
lattice spacing h. The contributions of higher-distance interactions are very small. The 
mutual separations and directions of the monopole currents corresponding to the couplings 
Qi are summarized in Table E We visualize the first seven most essential coupling constants 
in the monopole action in Figure |21 



coupling 


distance 


type 


coupling 


distance 


type 


91 


(0,0,0,0) 


^m(«) 


515 


(2,1,1,0) 


kfj,{s + 2fl + u + p) 


92 


(1,0,0,0) 


k^{s + jl) 


516 


(1,2,1,0) 


kfj.{s + p + 2i' + p) 


93 


(0,1,0,0) 


k^{s + z>) 


517 


(0,2,1,1) 


kfj,{s + 2i> + p + (t) 


94 


(1,1,0,0) 


kf^{s + fi + v) 


518 


(2,1,1,1) 


kfj,{s + 2fi + u + p + a) 


95 


(0,1,1,0) 


kfj.{s + z> + /5) 


519 


(1,2,1,1) 


kfj,{s + /i + 2z> + p + (t) 


96 


(2,0,0,0) 


k^{s + 2fi) 


520 


(2,2,0,0) 


kf,{s + 2p + 21)) 


97 


(0,2,0,0) 


k^{s + 2v) 


521 


(0,2,2,0) 


kf,{s + 2v + 2p) 


98 


(1,1,1,1) 


kn{s + fi + i' + p + a) 


522 


(3,0,0,0) 


k^{s + 3p) 


99 


(1,1,1,0) 


kfj_{s + fjL + V + p) 


523 


(0,3,0,0) 


k^{s + 3i>) 


910 


(0,1,1,1) 


k^{s + z> + p + (t) 


524 


(2,2,1,0) 


k^{s + 2jl + 21) + p) 


911 


(2,1,0,0) 


kn{s + 2/i + i>) 


525 


(1,2,2,0) 


k^{s + fi + 2i) + 2p) 


912 


(1,2,0,0) 


k^{s + jl + 2t') 


526 


(0,2,2,1) 


kf,{s + 2u + 2p + a) 


913 


(0,2,1,0) 


k^{s + 2i) + p) 


527 


(2,1,1,0) 


kp{s + 2p + 2i> + p) 


514 


(2,1,0,0) 


k^{s + 2fi + u) 









TABLE I: The quadratic interactions in the monopole action determined numerically. 



The action determined above takes into account all monopole trajectories. However, a 
typical monopole configuration in the confinement phase consists of one large monopole 
trajectory (percolating cluster) supplemented by a lot of small (ultraviolet) monopole clus- 
ters |j23|J. The percolating cluster fills the whole volume of the lattice and it makes a domi- 
nant contribution to the tension of the chromoelectric string. The properties of the largest 
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FIG. 2: The graphic representation of the first seven types of the quadratic interactions in the 
lattice monopole action schematized in Tabled 



percolating cluster were studied both numerically jl^, Hi 2m ^^'^ analytically 12; 
percolating cluster is associated with the monopole condensate j^, 1^ 
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FIG. 3: The couplings 171 and §2 of the monopole action (a) for all clusters and (b) for the percolating 
cluster. In this and other figures the error bars are smaller than the size of symbols and the scale 
b is shown in units of the string tension. 



If our determination of the monopole action is self-consistent then at large scales h the 
ultraviolet clusters should not give any contribution neither to the monopole action nor to 
the monopole condensate. The correctness of the first statement for the leading parameter, 
(?!, was confirmed in Ref. In Figures Efa), (b) we show the couplings gi and g2 for all 
clusters and for the percolating cluster. These couplings show an approximate scaling: they 
depend only on the product b = a-n and do not depend on the variables a and n separately 
when n > 3 are considered. The larger b the better scaling is. 

The comparison of the couplings computed on all clusters and on the percolating cluster 
only are shown in Figures |^a) and (b). Again, one can clearly observe that at large scales b 
the coupling constants evaluated on the different types of the monopole ensembles coincide 
with each other contrary to the small-6 case. 
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FIG. 4: The comparison of the couphngs (a) 51^2 and (b) 53^4 computed for all clusters and for the 
percolating cluster. 



V. MONOPOLE CONDENSATE FROM MONOPOLE ACTION 

To get the value of the monopole condensate we have to compare the monopole action 
calculated analytically in Section IIIII with the numerical results described in Section IIVI 
To this end we first note that due to the closeness of the monopole currents K^^^^ only the 
transverse part of the monopole operator (j28j) has a sense. Indeed, the shift of the quadratic 
operator S —>■ S + add' (with a being an arbitrary parameter) does not change the monopole 
action ()2fij) due to conservation condition (jZj). Therefore, in order to relate the theoretical 
and numerical results we need to get the transverse part of the operator ()28|). 

A simplest and, on the other hand, a practical way to do extract the transverse part 
of the quadratic monopole operator is to calculate the monopole action on a set of closed 
monopole trajectories, K^'^K We consider six types of such monopole trajectories which are 
depicted in Figure El 




FIG. 5: Set of lattice currents used to get transverse elements fi of the monopole action operator. 

Let us consider the analytical prediction for the transverse part of the monopole action. 
Since we are working in the fJ^ ^ q limit, we disregard 0{n~^) corrections to the quadratic 
action (j^H|) . The validity of such approximation is discussed below. The leading contribution 



12 



to the monopole action evaluated on closed trajectories K^^^ is 
where \K^'^^\ is the length of the trajectory K^'^'> and 



(36) 



do 


= (0,0,0) 






= 0.248028 


di 


= r'^^') (0,0,0) 


-P(3)(i^0,0) 




= 0.166665 


d2 


= (0,0,0) 


- (1^0,0) - 


ip(3)(2,o,0) 


= 0.181055 


ds 


= p(3)(o,o,0) 


- (1,0,0) - 


ip(3)(3,o,0) 


= 0.181292 


di 


= v(^\o,o,o) 


-p(3)(2,o,0) 




= 0.209836 


d5 




- (1,0,0) - 


|2^('Hl,l,0) 


= 0.176956 



(37) 



are the linear combinations of the values of the inverse three-dimensional Laplacian T)^^^ at 
certain points. The numerical values shown in Eq. (jHTj) correspond to the lattice 48'^. Below 
we call the combinations /'•*•' of the couplings as "transverse couplings". 

Using Table m one can get the transverse combinations of couplings corresponding to the 
numerically calculated action: 



/o = gi + 2g2 + 2gG + 2g22 

/2 = gi + 1(^2-^3-^4) - 

U = gi + g2- gr- gu , 



fi = gi- gs, 

h = gi + g2- 1^3 - gA 

fb = gi- Igs - lg5 ■ 



l96 



\9ii 



\g23 



(3J 



Note that the transverse components of the analytical action with two free param- 
eters should describe six transverse combinations ()38|1 obtained numerically. We fit the fi 
components by ()36|1 independently for each i = 0, 1, . . . , 5 and then compare in Table |n] 
the fitting parameters rj and muv as a self-consistency test. Since we are working in the 
fi ^ 1 limit we fitted the numerical data for the n = Q blocked monopoles. A lower value 
of n corresponds to the smaller scale b and in this case we notice sizable deviations of the 
numerical results from our fitting function. This is expected because we are working in the 
limit b oo. One the other hand, the higher value, n = 8, correspond to the small lattice 
size of the coarse lattice, (N/n)'^ = 6^, which may lead to large finite-volume artifacts. 
Therefore we concentrated on = 6 blocked monopoles. 

The fits of the transverse couplings of the monopole action corresponding both to all 
monopole clusters and to the percolating cluster are visualized in Figures Efa) and (b), 
respectively. The best fit parameters obtained from the fits of different transverse couplings 
fi (Table |n} are very close to each other what provides a nice self-consistency test of our 
approach. Moreover, the value of the monopole condensate t] calculated in large-6 limit 
from the all-clusters and the percolating cluster monopole action are the same within error 
bars, as expected. The numerical value of the monopole condensate (obtained by averaging 
of the results of the six independent fits) is ?7 = 243(6) MeV. 

Finally, let us discuss the validity of the large /i approximation used in this paper. We 
are working in the range of momenta by/a ~ 1 ... 4. The mass of the dual gauge boson 
obtained from the fitting of the string profile by a classical string solution is estimated 
as Mb ~ 1.3 GeV ~ 3^/a. Therefore the value of fi, Eq. (j2SI), is in the range ~ 3 ... 12. 
There are two types of corrections to our analytical results: (i) the exponentially suppressed 
corrections to the operator J^~^ (discussed in Appendix (jB))) are smaller than 5%; (ii) the 
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coupling 








all clusters 


max cluster 


all clusters 


max cluster 


fo 


0.521(25) 


0.509(23) 


0.046(9) 


0.042(8) 


h 


0.577(41) 


0.580(45) 


0.020(9) 


0.022(10) 


J2 


0.565(34) 


0.537(37) 


0.031(9) 


0.025(9) 


h 


0.544(32) 


0.522(35) 


0.032(9) 


0.026(9) 


h 


0.554(28) 


0.532(30) 


0.041(9) 


0.034(9) 


h 


0.591(38) 


0.590(42) 


0.025(9) 


0.026(10) 


average 


0.552(13) 


0.534(13) 


0.036(4) 


0.031(4) 



TABLE II: The values of the condensate r/ and the ultraviolet cutoff tuv obtained in a set of 
independent fits of the n = 6 transverse monopole couplings H38() by function (|36|) for the all 
monopole clusters and for the percolating monopole cluster. The best parameters of the overall fit 
of /i . . . /g are shown in the last row. 




FIG. 6: The fits of the transverse n = 6 monopole couplings (|38j) by function (|36j) (a) for the all 
monopole clusters and (b) for the percolating monopole cluster. 



0(/i ^) correction of Eq. is of the order of 10% because of the local nature of the Q 
and Ti operators, and due to low values of the inverse Laplacian, V^^\Q,Q,Q) ^ 1/4. Thus 
we estimate the systematic corrections to the value of the monopole condensate to be of the 
order of 15%. Taking into account the systematic errors we get r] = 243(42) MeV. 



VI. DISCUSSION AND CONCLUSION 



We have obtained the value of the monopole condensate using the method of blocking 
from continuum to the lattice. Namely, we have obtained numerically the effective monopole 
action in the Maximal Abelian projection of quenched SU(2) lattice QCD. Then we have 
calculated analytically the effective lattice monopole action starting from the continuum dual 
Ginzburg-Landau model. In our simulations we restricted ourselves to the large values of the 
parameter b. This parameter defines a scale at which the monopole charge is measured on the 
lattice. In large-6 limit the action of the monopoles is dominated by the quadratic part, and 
higher monopole interactions are suppressed. Thus in our analytical calculations we have 
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neglected the quantum contributions of the scalar monopole fields which are responsible for 
the higher order corrections to the effective monopole lagrangian jl8l |. 

The comparison of the numerical and analytical results for the blocked action gives us 
the value of the monopole condensate, t] = 243(42) MeV. This value is in a quantitative 
agreement with another estimation of the monopole condensate, t] = 194(19) MeV, obtained 
in Ref. using a completely different method. Moreover, we have shown that our method 
is self-consistent, since is allows to describe various quadratic interaction of the monopole 
action using approximately the same values of the monopole condensate. 

A few words about the ultraviolet cutoff tuv are now in order. Thus cutoff - which enters 
the effective monopole action ()28|) - is an independent fitting parameter of the effective 
monopole action at large scales, Eq. (jHUj) . In this paper we have neglected the fluctuations 
of the monopole scalar flelds since we were working at large scales b. Effectively, this cor- 
responds to taking the London limit of the Ginzburg-Landau model. The London limit 
possesses known logarithmic divergences {i.e., the tension of the Abrikosov vortex is log- 
arithmically divergent function of an ultraviolet scale). The physics of the monopole fleld 
fluctuations is "hidden" in the value of this cutoff. Strictly speaking, we have to renormalize 
the model and consider the monopole fleld fluctuations to relate a logarithmic divergence to 
the values of the physical parameters entering the lagrangian of the model. This procedure 
becomes meaningful at small scales b. 

At small values of the scale b the higher order interactions ( four-point, six-point, etc) 
become essential Thus at short distances the scalar monopole fleld contribute to the 
effective monopole action. From the point of view of the blocking from continuum, at small 
values of b the couplings of the monopole action become dependent on the parameters of the 
potential of the monopole fleld. Thus, a comparison of the effective monopole action with 
the blocked action at small scales b may allow us to determine the shape of the monopole 
potential. We discuss this problem in a forthcoming publication |i28i] . 
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APPENDIX A: PROOF OF CLOSENESS OF LATTICE MONOPOLE CURRENTS 

In order to prove the relation (|7j) it is convenient to represent the lattice monopole 
current (0) as the integral over momentum. Using Eq. ()20|) and Eq. (j^^ we get: 

^ ^ap{C^,s, x) = i{p^, S^^ - 6^^) b^ Q^{pb) e'^^'^P''^ , (Al) 

where Sq/3(C, x) = S^^(x), the vector Qa is given in Eq. (fT^ . and no summation over index 
7 is assumed. Then Eq. (jS} can be rewritten as follows: 

Ks,, = {P, 5,, - Pu 5,,) ki-p) ^ QM) ^'^'^'''^ 
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= ki-P) Q.ipb) e-''^^''^ , (A2) 

where k^{p) = J dx k^{x) e~''^^'^^ is the Fourier transformed continuum monopole current. 
There is no summation over index 7 in Eq. ()A2|) . To get the second hne of Eq. ()A2|) we used 
the closeness condition of the continuum monopole currents, 

pMp)=^- (A3) 

According to Eq. (jH) the lattice monopole currents, i^s,^, are associated with the centers 
of the three-dimensional cubes C^^^. The positions of the cube centers are characterized 
by the integer-valued coordinates s. The corners of the cubes belong to the original lattice 
while the monopole currents themselves are associated with the dual lattice. The sites of 
the dual lattice are shifted by the four-dimensional vector (1/2, 1/2, 1/2, 1/2) with respect 
to the sites of the original lattice. Thus, the center of the cube Kg^^ does not belong to the 
dual lattice because the coordinate of the center of the cube corresponds to the time-slice 
of the original lattice. In our coordinates, the monopole current defined on the cube i^^^^ 
must be associated with the point *s = s + fl/2 belonging to the dual lattice. 

Thus, the closeness condition ((7j) at the site *s of the dual lattice is 

4 4 

7=1 7=1 
= 2ib' J ^ ~k,{-p) sm{bp,/2) Q,{pb) e-^^(P'^) . (A4) 

Using Eq. (fTIJ|) we notice that 2sm{bp^/2) Q^{pb) = p-^Qipb), where the quantity Q{x) = 
^^x )-! does not carry any Lorentz index. Then Eq. ()A4|) together with the conservation of 
the continuum monopole charge ()A3|) implies the closeness of the lattice monopole currents, 

{d'K), = lb' I -^^p;k,{-p) Q{pb) e-^"^^''^ ^ . 



APPENDIX B: CALCULATION OF THE OPERATOR T'^ 

In this Appendix we calculate the expression for the inverse operator JF^^, presented in 
Eq. (121, for /i = MBb > 1. 

Let us consider first the diagonal components of the inverse operator JF~^. Without loss 
of generality we take = z/ = 4 and s' = 0. We get: 

J (27r)4 qi + q^ + jj,^ fj^ \ qi/2 J 

It is convenient to introduce the additional integral, 

1 



/I + + /i^ 



POO 

/ dte-(«*+^^+'^')*, 
Jo 



(B2) 
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and represent the integral ()B1|) in the form: 



°? 3 3 

^:oM = / dte-'^'*Po(54,t) J]Pi(s„t)n^2(s,,t) 



(B3) 



3 = 1 



where 



Po(s,t) 
Pi(s,t) 



dg 



1 



(B4) 



-oo 
oo 



27r 2Vvrt 



sm - ■ e 
27r 2 



2V7rt 



-oo 
oo 



dg /sing/2\' 
2^ V g/2 J 



1 

+ 2 



+ {s 



l)Erff 



s - 1 
V 2v^ 



- 2sErf 



2v^ 



(B6) 



and Erf(x) is the error function, 



Erf(x) = ^ / dye'J' . 



To calculate the off-diagonal components of the inverse operator JF ^ we take (yU, v) 
[1, 2) and s' = (again, without any loss of generality): 



-^50,12 



d^^g -gig2 TTsingj/2 T-r/singj/2 



n 



(2ir)-' ,2 + ^2 11 ,_/2 i-1 V ,j/2 

/•OO 2 4 

/ dte-'^'*n^3(s.,t)-n^i(^i'^) 

i=l ,=3 



(B7) 



where 



P3(s,t) = 2 



dq . g 



sm - ■ e 

27r 2 



-tq +iqs 



2V7rt V 



/g-(2s-l)Vl6t _ g-{2s+l)Vl6t 



[B8) 



Eqs. ()BlllB8|) represent the exact expressions for the diagonal and off-diagonal elements 
of the inverse operator JF~^. Unfortunately, due to the presence of the Erf -functions in 
Pi, Eq. ()B5|) . the integrals ()B1|) and ()B7|) can not be taken analytically. However, in the 
limit /i — > oo, which corresponds according to Eq. ()25|) to large blocking scales 6, leading 
contributions to these integrals can be easily estimated. 

Let us first consider Eq. ()B1|) . The main contribution to this integral is coming from the 
region of small t. At small t the error function can be represented as 



Erf(x) = sign(x) (1 + Ofx"^)) for a; S> 1 

^J'KX 



(B9) 
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Therefore at general values of s the expression ()B1|1 is given by a sum of integrals of the 
form 

POO 

I{fx,s)= dt exp{-fxH- 3'^/t + C \ogt} , (BIO) 
Jo 

where C is a constant of the order of unity and the quantity s depends on the value of s 
[i.e., s = s/2, {s — l)/2, etc.). The value of s is either of the order of unity or zero. At 
§ ~ 1 and large /i we get J(/i, s) ~ exp{— 2yus} -C 1. Thus the integral ()B10|) with s 7^ is 
exponentially suppressed and therefore it will be neglected below. The leading contribution 
to the operator J-'~^ comes from the integrals of the form ()B10|) with s = which are 
saturated at small t. 

Using the expansion ()B9j) we get in the leading order in the limit t — > 0: 

Pois,t) = -^5, + 0(e-=-^*/*), (Bll) 

2v TTt 



Pi{s, t) = ^ A, + 0(e-=°°^*/*) , (B12) 

2v TCt 



P2{s, t) = ^^As + 6s + 0(e-^°'^^*/*) , (B13) 
P3{s,t) = 0(6-^°=^^*/*), (B14) 



where 



1, s = 
0, otherwise 



1, s = l,-l 
As={ -2, s = , (B15) 
0, otherwise 



are the Kronecker symbol and the one-dimensional lattice Laplacian, respectively. 

According to Eqs. ()B7IB14|) . the elements with fi ^ u of the operator JF~^ are exponen- 
tially suppressed, ^^^^ ~ 0{e~^°'^^^'^). As for the diagonal elements of this operator, /i = i^, 
we get 

00 33 

J'^OM = / E n (^«. + vf'^^O + 0(e-=°°^*'^) (B16) 

^ { 1=1 j=i ^ 



1 e 



r(0,tc/y/i^) (^As.ds^Ss^ + cyclic) + - (^A^^A.^^.g + cyclic) 

0/ -constMN _ 



Here T is the incomplete gamma function, r(a, x) = t"" ^ e * dt, and "cyclic" means cyclic 
permutations over the indices st. To get Eq. ^BT^ we used Eqs. 1IB12IK18jl . We 

also introduced the ultraviolet cutoff, t[/y, to regularize the logarithmically divergent piece 
of Eq. ()B16|) . Noticing that V'^\s} = As^6s2Ss3 + cyclic is the three dimensional Laplacian 
we get the final expression for JF~^ presented in Eq. fl27p . 
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